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Abstract 

A method for determining permeability of phospholipid bilayer based on the 
osmotic swelling of micrometer-sized giant unilamellar vesicles (GUVs) is 
presented as an alternative to the two established techniques, dynamic light 
scattering on liposome suspension, and electrical measurements on planar 
lipid bilayers. In the described technique, an individual GUV is transferred 
using a micropipette from a sucrose/glucose solution into an isomolar solution 
containing the solute under investigation. Throughout the experiment, vesi- 
cle cross-section is monitored and recorded using a digital camera mounted 
on a phase-contrast microscope. Using a least-squares procedure for circle 
fitting, vesicle radius R is computed from the recorded images of vesicle 
cross-section. Two methods for determining membrane permeability from 
the obtained R(t) dependence are described: the first one uses the slope of 
R(t) for a spherical GUV, and the second one the R(t) dependence around 
the transition point at which a flaccid vesicle transforms into a spherical one. 
We demonstrate that both methods give consistent estimates for membrane 
permeability. 
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1. Introduction 

The cell membrane physically separates the cell interior from the envi- 
ronment. The membrane is selectively permeable: it permits the passage of 
some species of molecules or ions to- and from the cell, while it blocks the 
transport for other species. Around 1900, while studying the working of gen- 
eral anaesthesia, Hans Meyer in Marburg [lj and Ernest Overton in Zurich 
[2| independently devised a simple rule to predict membrane permeability. 
They established that the more a molecule species is soluble in lipid, the 
greater the cell permeability for this molecule species is. While this rule can- 
not account for transport processes not known at that time, such as those 
mediated by membrane carriers, channels (e.g., aquaporins (3)]) or pumps, 
nor does it acknowledge membrane inhomogeneities, such as rafts Q present 
in the biological membrane, the Meyer-Overton rule seems to withstand the 
test of time [HI, . At the same time, there is a continuing interest in both 
basic physics of general anaesthesia [7| and passive membrane permeability 

While a significant corpus of publications on permeability properties of 
biological cells has been accumulated since Meyer and Overton, experimental 
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access to an isolated lipid bilayer only became available in the early 1960s, 
when a technique for preparation of thin films separating two aqueous com- 
partments has been devised, known as the black lipid membrane [ij]. An 
arguably even more influential technique which was developed at roughly 
the same time allowed for producing artificial lipid vesicles, or liposomes 



10(. Vesicles are osmotically sensitive structures which swell and shrink in 
response to changed osmotic conditions. Even though the sub-micrometer li- 
posomes which are the easiest to produce cannot be visualized directly, their 
size can be estimated via light scattering (see [l2[ for a review on the 
optical methods in determining membrane permeability). 

Both model lipid bilayer systems - planar lipid bilayers and liposomes 
have proved extremely fruitful in the studies of membrane permeability 



13J. Planar lipid bilayers are well suited to electrical characterization, as 
the two chambers separated by the bilayer are both easily accessible, which 
allows simple placement of macroscopic electrodes. Finkelstein 14( mea- 
sured the permeability of a planar lipid bilayer for water and seven other 
non-electrolytes in an attempt to resolve the mechanism by which neutral 
molecules and ions permeate the membrane. One possibility is the solubility- 
diffusion mechanism, which assumes that the permeating species dissolves in 
the hydrophobic membrane, diffuses across, and leaves by redissolving into 
the other aqueous compartment. Another possibility is that permeation oc- 
curs through hydrated transient defects, which appear as a result of thermal 
fluctuations. That study concluded that both water and non-electrolytes 
cross the membrane through the solubility-diffusion mechanism. A similar 
but more comprehensive study was repeated a few years later by Orbach 



and Finkelstein [15(. Walter and Gutknecht [16[ examined the correlation 
between the membrane permeability for 22 solutes and their partition coeffi- 
cient between water and any of the four examined organic solvents, and found 
a very high correlation with hexadecane and olive oil, and a less pronounced 
correlation in the case of octanol and ether. 



In an early work [17], Bangham and coworkers made use of the fact that 
the total volume of liposomes in the suspension is proportional to the re- 
ciprocal of the optical extinction, which allowed them to examine the per- 
meability of membrane for water and various solutes via osmotic swelling 



and shrinking of liposomes. As early as 1933 Jacobs noticed [18| that the 
volume of a cell transferred into a solution of permeant solute transiently de- 
creases, reaches some minimal value, and then starts increasing. Sha'afi and 
coworkers [l9[ employed this phenomenon and the Kedem-Katchalsky for- 
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malism 20] to compute the permeability of erythrocyte membrane for urea. 
De Gier and coworkers used the initial slope of the reciprocal of the optical 
extinction for determining membrane permeability of liposomes 211. Hill 



and Cohen 22j brought the "minimal volume" technique to the experiment 
with liposomes as well. A comprehensive review of the use of liposomes in 
membrane permeability studies is given by de Gier 23]. Using established 
techniques, Paula et al. [24| did an extensive study in another attempt to 



resolve the standing debate between the solubility-diffusion mechanism and 
the hydrated transient defects as a primary pathway, and concluded that 
except for ion permeability of lipid bilayers composed of phospholipids with 
short chain lengths, solubility-diffusion mechanism seems to be the domi- 
nant effect. Examining the known phenomenon that upon transfer into a 
hypotonic medium, vesicles swell, and, if the gradient is large enough, burst 
and expel part of their content, Shoemaker and Vanderlick studied the influ- 
ence of membrane composition on the extent of leakage 25], and found out 



a correlation between the membrane resistance to burst and its stretching 
modulus. 

Both planar lipid bilayers and liposomes as model bilayer systems have 
their drawbacks, too. The original "brush" technique of producing planar 
lipid bilayers has been limited by pockets of residual solvent trapped be- 
tween the two bilayer leaflets, which affects membrane properties. While an 
improved deposition method 26| virtually eliminated this problem, the lim- 
ited lifetime of the membranes- most often less than one hour - remains a 
persisting problem which limits the duration of the experiment. Liposomes, 
on the other hand - in particular large unilamellar vesicles (LUVs) with a 
diameter of 100-200 nm - have proved to be extremely stable. The interpre- 
tation of dynamic light scattering (DLS) experiments of osmotic shrinking of 
LUVs is not trivial, as it involves the transformation of shapes predicted by 
the area difference elasti city model (27| to the hydrodynamic radius, which 



is characterized by DLS [28|, |29|. Also, it has been argued that LUVs, which 
consist exclusively of high-curvature membrane regions, serve as a poor model 
of biological cell membranes. An attempt to resolve the possible dependence 



of permeability on the membrane curvature [30J was inconclusive, as the au- 
thors ascribe the observed differences in permeability to the problems they 
experienced with planar lipid membranes. Finally, an effect which occurs 
in both systems, but is more prominent with the planar lipid bilayer, is the 
unstirred layer effect 311 ] . In general, the concentration of solute adjacent 
to the membrane differs from its concentration in bulk. It is the concen- 
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tration of solute immediately adjacent to the membrane which determines 
the permeation of solute across the membrane, while the concentration in 
bulk is the one that is usually known. In both the experiment design and 
the interpretation of the experimental findings, one needs to be aware of this 
discrepancy. 

While both planar lipid bilayers and liposomes have been used as model 
membrane systems for studies of membrane permeability since 1960s, stud- 
ies employing GUVs, which allow for a direct visualization of the process, 
appeared decades later [32|, chiefly due to a lack of suitable techniques for 
preparation and manipulation of GUVs in those early days. In the paper 
by Boroske et al. 32J, the authors describe the experiment in which GUVs 
were prepared in pure water and subsequently transferred into a solution of 
either glucose or NaCl (concentrations used ranged from 1.5-20 mM) while 
their size was monitored using phase-contrast microscopy. Upon transfer, 
vesicles shrank in size; the process of shrinking depended on the vesicle size. 
Vesicles with radius R < 10 (J.m shrank with a linear time dependence: 
R(t) = Ro — V w PAct, where V w is the molar volume of water, P is membrane 
permeability for water (water filtration coefficient), Ac the solute concentra- 
tion difference, and t time. From measured data, the authors inferred the 
water filtration coefficient, P = 41 |a.m/s. Larger vesicles (r > 10 \xm) under- 
went a phase of "instability" in which the vesicle was flaccid, and after which 
a spherical shape was reestablished. The authors dismissed the idea of dis- 
solving lipid molecules into the outer medium as a plausible explanation for 
the apparent decrease of vesicle surface area, in particular since they also no- 
ticed formation of smaller satellite spherical vesicles, seemingly connected to 
the mother vesicle. Instead, they proposed a mechanism of concerted flipping 
of lipid molecules from the inner membrane leaflet into the outer membrane 
leaflet, induced by the flow of water. Recently, membrane permeability has 
also been studied of GUVs made of block copolymers [33|, [34 . 



The remaining of this chapter is structured as follows. First, the exper- 
imental section introduces the system and presents the immediate experi- 
mental results. A section on the theory of membrane permeability offers a 
review of the few selected phenomenological models for membrane permeabil- 
ity, with a special emphasis on the influence of membrane elasticity and the 
swelling-burst cycle. The section concludes with a theory of the continous 
transition between the "ironing" and the stretching regimes of the osmotic 
swelling of a vesicle. The section on experimental analysis demonstrates the 
calculation of membrane permeability based on the swelling-burst cycle, and 
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compares its result with the calculation based on the analysis of the transi- 
tion between the "ironing" and the stretching regimes. We conclude with a 
discussion of the merits and limitations of the presented method. 



2. Experimental section 

2.1. Materials and methods 

D-(+)-glucose, D-(+)-sucrose, glycerol, urea, and ethylene glycol were 
purchased from Fluka (Buchs, Switzerland). Methanol and chloroform were 
purchased from Kemika (Zagreb, Croatia), l-palmitoyl-2-oleoyl-sn-glycero- 
3-phosphocholine (POPC) was purchased from Avanti Polar Lipids (Alabaster 
USA). All the solutions were prepared in double-distilled sterile water. 

A suspension of POPC GUVs in 0.1 or 0.2 mol/L 1:1 sucrose/glucose so 



lution was prepared usin g an electroformation method, described in Ref. |35 



with some modifications [36|, [37| • Lipids were dissolved in a mixture of chlo- 
roform/methanol (2:1, v/v) to a concentration of 1 mg/mL. A volume of 
25 \xL of the lipid solution was spread onto a pair of Pt electrodes and dried 
under reduced pressure (water aspirator; « 60 mmHg) for 2 hours. The elec- 
trodes were then placed into an electroformation chamber, which was filled 
with 0.1 or 0.2 mol/L sucrose. AC current (8 V, 10 Hz) was applied, and the 
voltage and frequency were reduced in steps to the final values of 1 V and 



1 Hz [37| • Subsequently, the chamber was first drained into a beaker and then 
flushed with an equal volume of isomolar glucose solution, thus resulting in 
a suspension of GUVs containing entrapped sucrose in a 1:1 sucrose/glucose 
solution, which increases the contrast in a phase contrast setup and facilitates 



vesicle manipulation [38J. This procedure yields mostly spherical unilamellar 
vesicles, with diameters of up to 100 \mv. 

An inverted optical microscope (Nikon Diaphot 200, objective 20/0.40 
Ph2 DL) with micro-manipulating equipment (Narishige MMN-l/MMO-202) 
and a cooled CCD camera (Hamamatsu ORCA-ER; C4742-95-12ERG), con- 
nected via an IEEE- 1394 interface to a PC running Hamamatsu Wasabi 
software, was used to obtain phase contrast micrographs. In the streaming 
mode, the camera provides 1344 x 1024 12-bit grayscale images at a rate of 
8.9 images/s. 

In the experiment, an individual spherical GUV is selected, fully aspirated 
into a glass micropipette whose inner diameter exceeds the vesicle's diameter, 
and transferred from a solution containing solutes of very low membrane 
permeability (1:1 glucose/sucrose) into an iso-osmolar solution of a more 
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Figure 1: A sequence of micrographs showing an individual vesicle burst. Vesicle radius 
before the burst was 32.69 urn; immediately after the burst, it shrunk to 30.47 urn. The 
images were taken at 0.11 s intervals. The bar represents 20 urn. 



permeable solute (glycerol, urea, or ethylene glycol), where the content of 
the micropipette is released, and the micropipette is subsequently removed. 
Vesicle response is recorded using a CCD camera mounted on the microscope. 

2.2. Experimental results 

Upon transfer into a solution of permeating solute, vesicles start to swell 
until the membrane critical strain is reached. At that point, the membrane 
ruptures, and the vesicle ejects part of its internal solution (figured]). Af- 
ter the burst, the membrane reseals, and another cycle of swelling starts. 
The observed sequences of swelling-burst cycles ranged from 3 to over 40 
successive bursts, with an average in our sample being 15.4. 

For a quantitative analysis, the radius of the vesicle cross-section in each 
recorded micrograph was determined with a GNU Octave script using a least- 
squares procedure for circle fitting [39j . Figure [2] shows a typical time course 
of vesicle radius after transfer. In about 80% of all cases, we observed a 
transient maximum of radius shortly after the transfer. We attribute this 
phenomenon to a slight hypertonicity of the target solution. This causes the 
vesicle to deflate rapidly, the rate of volume change being determined by 
the vesicle size and the permeability of phospholipid membrane for water. 
A deflated vesicle changes its shape from spherical into a shape which can 
be approximated with an oblate spheroid. Small deviations from this shape, 



which originate from the effects of gravity |40j, are neglected here. Due to 
gravity, the axis of rotational symmetry is aligned with the vertical, which, 
in our experimental setup, also coincides with the direction of the optical 
axis. Thus, the observed cross-section radius is the oblate spheroid longer 
semiaxis, and consequently it increases with the decreasing volume, while the 
membrane area remains unchanged. 
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Figure 2: Time dependence of the radius of a vesicle cross-section upon transfer from a 
0.1 mol/L 1:1 glucose/sucrose solution into a 0.1 mol/L solution of ethylene glycol. Clearly 
visible are the initial maximum, followed by a saw-tooth pattern of repeated swelling-burst 
cycles. 



Concurrently with the initial vesicle deflation due to the slight hyper- 
tonicity, diffusion of permeating solute into the vesicle interior, accompanied 
by the osmotic influx of water, also takes place. The rate of this latter process 
is however dominated by the permeability of the membrane for the given so- 
lute (in our case, glycerol, urea, or ethylene glycol), and is thus much slower 
than the initial rapid deflation due to the mismatch in tonicity. With time, 
however, the initial trend of vesicle deflation is overturned and the vesicle 
starts inflating again. Inflation can be qualitatively divided into two phases. 
In the first phase, the vesicle is still partially deflated, and the influx of solute 
and the accompanying osmotic influx of water increases the vesicle volume, 
all while the membrane area remains unchanged. This "ironing" phase is 
characterized by the decrease of the radius of vesicle cross-section. At some 
point, the vesicle reaches a spherical shape, which can be observed as a local 
minimum of the cross-section radius. 

Further permeating of solute into the vesicle interior, accompanied by 
the osmotic influx of water, causes the vesicle to inflate while maintaining 
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a spherical shape. In this phase, the increase of vesicle volume can be ob- 
served as an increase of vesicle cross-section. In order to accomodate the 
increased vesicle volume, the membrane needs to stretch. This stretching 
process continues until the critical strain for the membrane is reached. At 
this point, the membrane ruptures, and part of the vesicle interior is ejected 
outside. When the vesicle volume is thus reduced, the membrane reseals 
again. As the concentration difference for the permeating solute persists, 
which serves as a driving force for the diffusion of the permeating solute into 
the vesicle interior, this means that at the same moment, the vesicle volume 
starts increasing again, and another swelling cycle commences. Repeated 
swelling-burst cycles give yield to the characteristic saw-tooth pattern, when 
the vesicle radius is plotted against the elapsed time (figure [2]). 

In total, 47 recordings of vesicle transfer from a 1:1 sucrose/glucose so- 
lution into an isomolar solution of glycerol, urea, or ethylene glycol were 
selected for further analysis. Out of these, 15 transfers were into glycerol 
(5 at 0.1 M and 10 at 0.2 M), 15 transfers into urea (8 at 0.1 M and 7 at 
0.2 M), and 17 transfers into ethylene glycol (12 at 0.1 M and 5 at 0.2 M). 
Figure [3] shows the critical strain (v4 crit — A )/A = (R crit /R ) 2 — 1, where 
Ro is the radius of a relaxed spherical vesicle, and R cr it is the radius of a 
critically strained vesicle. The median critical strain obtained for a total of 
738 recorded vesicle bursts is 0.033, and mean critical strain is 0.038 with 
standard deviation 0.024. 



A published value for the critical strain is around 0.04 [41], |42j. Our own 



estimate based on the first burst in the sequence alone is a little higher: 



0.055 ± 0.02 [43]. This may indicate that the membrane might not always 
perfectly reseal, and that a local defect present in the membrane makes it 
more likely to rupture at a lower strain. 

3. Theoretical section 

In this section, we will briefly review a few theoretical models used for 
the description of solute permeation across cell membrane, starting from the 
most general Kedem-Katchalsky formalism, based on nonequilibrium ther- 



modynamics [44] . then showing how in the case of separate pathways for 
water and solute transport it can be simplified into a two-parameter model, 
and subsequently when a two-parameter model can be simplified into a one- 
parameter model in the case when the membrane permeability for water 
greatly exceeds its permeability for solute. We will discuss the influence of 
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Figure 3: A histogram of experimentally determined critical strain, defined as (R C rit/Ro) 2 — 
f, where Rq is the initial and i? cr it the final, or critical, radius for each of the swelling 
cycles. 



finite membrane elasticity on the apparent membrane permeability, review 
the work on the repetitive swelling-burst cycles, and conclude with the theory 
of the continuous transition between the "ironing" and the stretching regime 
of a vesicle which swells osmotically due to the permeation of a permeable 
solute. 



3. 1 . Kedem-Katchalsky formalism 

Starting from the principles of non-equilibrium thermodynamics, Kedem 



and Katchalsky arrived at a model with three parameters 20]. A simple 
qualitative argument can be offered why no less than three parameters are 
required for the description of the transport of solvent and solute across 
a membrane. In the simpler case of free diffusion of a single solute in a 
solvent, molecules of solute move relative to the molecules of solvent, and a 
single parameter - diffusion coefficient - is required to describe the friction 
between the solvent and the solute. In the case of membrane transport, two 
additional coefficients are required, one describing the friction between the 
molecules of solvent and the membrane, and another describing the friction 
between the molecules of solute and the membrane. 

As their starting point, Kedem and Katchalsky [20| choose the rate of 
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entropy production during the permeation of solute and water across the 
membrane separating the interior compartment from the exterior, 



Here, /i s and fi w denote the chemical potential of solute and solvent (water), 
respectively, and diV*/dt and dN^/dt denote the number of moles of solvent 
and water entering the interior compartment per unit time. 

Formulating Eq. ([T|) in terms of a dissipation function per unit area, 
$ = (T/A) dS/dt, one obtains 

$ = fat, - /4) < + (/^ - /4) < • (2) 

Here, we have denoted n w = (1/A) dN^/dt for water and accordingly for the 
solute. Eq. (j2J) is a particular case of dissipation function, which in general 
assumes the form $ = JjXj, where Jj represent flows and X { represents 
the corresponding conjugate "forces". In Eq. (j2J), the differences in chemical 
potential act as conjugate forces. We shall now rewrite Eq. (j2J) in a way 
where the driving forces will be quantities easier to evaluate experimentally. 
For an ideal solution, the difference in chemical potential can be expressed 

as 

fj, e -// = V Ap + aAlnx , (3) 

where V denotes the partial molar volume of the constituent, x its mole 
fraction, Ap the pressure difference between the external and the internal 
compartment, and a = RT, with R being the molar gas constant and T 
the absolute temperature. Assuming a dilute solution, where the volume 
fraction ip of the solute is low, tp = c s V s 1, with c s = N s /V w being the 
molal concentration of solute, one can rewrite Eq. d3J) for the solute: 

Ac 

^-S s = V s Ap + a^ , (4) 

C s 

where Ac s = c e s — d s is the difference and c s = (c e s + c*)/2 is the average of the 
concentration of solute in both compartments. An analogous relation can be 
written for the solvent, 

a4,-/4, = ?«Ap-«— > (5) 
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with c w = (1 - ifi)/V w « 1/K,. Introducing Eq. (0} and Eq. © into Eq. © 
and rearranging, one arrives at another expression for the dissipation func- 
tion: 

— - — )aAc s . (6) 

Cs C s J 

In Eq. ([6]), $ is expressed in terms of forces commonly used in the perme- 
ability studies: X v = Ap is the hydrostatic pressure, and Xp = a Ac s is 
the osmotic pressure. The conjugate flows are J v = h w V w + n s V s , the total 
volume flow per unit area, and Jd = n s /c s — h s /c s , the relative velocity of 
solute with respect to solvent, which serves as a measure of exchange flow. 

It is assumed that each flow J present in the system in general depends 
on all the forces X acting in the system, and if the forces are sufficiently 
small, the relationship is linear: 

J\ = h\\X\ + L12X2 

J2 = Lii\X\ + L22X2 



Here, are phenomenological transport coefficients. The Onsager reci- 



procity relations [45|, |46[ state that the matrix of transport coefficients is a 



diagonal one, L ki = L ik . 

Applying this formalism to our case, we obtain 

J v = L p Ap + L p d a Ac s , (7) 
J D = L pD Ap + L D a Ac s , (8) 

where we have already taken into account Ld p = L P £>. The second law of 
thermodynamics requires that the diagonal terms are non-negative, while the 
off-diagonal terms are only constrained by the relation L v Ld — l? pD > 0. It 
is the off-diagonal terms though which are responsible for a concentration 
difference producing a volume flow, or vice versa. 

For easier comparison with the experimental results, it is convenient to 
transform {L p , Ljj, L p £>} to another set of coefficients. One of them is the 
reflection coefficient a, introduced with 

L p d = -crL p . (9) 

Two special cases which can be considered include a non-selective membrane 
(cr = 0) and an ideally selective membrane, permeable only for the solute 
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(er = 1). Another parameter commonly denned is the mobility of the solute 
uj: 

L P L D -L D 

u = c s = (L D - L p o )c s . (10) 

Lp 

It can be shown that u is chosen in such a way that for permeability mea- 
surement at constant volume (J v = 0), one can write n s = uRT Ac s . 

Expressing the volume flow and the exchange flow in terms of {L p , u, a}, 
one obtains 

J v = LpAp — aLpRT Ac s , (11) 
h s — (1 — cr)LpC s Ap + [oj — er(l — a)L p c s ] RT Ac s . (12) 

Eq. ( ITTj) and Eq. (TT2"j) obtained above pertain to a two-component system, 
in which binary solutions of the same solvent and solute are separated by a 
membrane. In order to allow for a comparison with the experimental results, 
we need to consider a slightly more complicated system, comprising of a 
membrane separating two ternary solutions of the same solvent and two 
different solutes. The membrane is permeable for the solvent and one of the 
solutes, but impermeable for the other solute. As in the previous example, 
we treat the solution as ideal. 

For a dilute solution, the difference in the chemical potential for the 
solvent (Eq. flSJ)) can now be written as 

Afi w = —a— — - a^^- + V w Ap, (13) 

c w c w 

where the index s pertains to the permeating solute, and the index n to the 
non-permeating solute. While Eq. fll]) remains valid, we have to consider 
the contribution of the non-permeating solute to the difference in osmotic 
pressure, 

ATl n = aAc n . (14) 

Considering the same flows as before, one can write the expression for the 
dissipation function Eq. (jSj), and when substituting into it the expressions for 
the difference in the chemical potential for the solute, Eq. (Tjf]), and solvent, 
Eq. (fT3|) . one obtains the expressions for the forces conjugate to the flows: 

X v = Ap — AU n , (15) 
X D = a Ac s + Lp AH n . (16) 
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Using the expressions ( TT51fTt3|) . one can obtain expressions analogous to (jllll2p . 
obtained in the case of a binary solution: 

J v = L p {Ap - AU n ) - aL p (a Ac s + <p AU n ) , (17) 
n s = (l- a)L p c s (Ap - ATl n ) + [u- <r(l - a)L p c s \ (a Ac, + ip ATl n ) . (18) 

A more condensed expression of Eq. ({TBI) can be obtained if Eq. (|T7|) is taken 
into account: 

h s = (1 — cr)c s J v + cu(a Ac s + ip AIi n ) . (19) 

In dilute solutions, <p AII n is often negligible with respect to a Ac s . On the 
other hand, the contribution AII„ is important. In many biologically relevant 
experiments, Ap = while AII n ^ 0. 

3.2. Two-parameter model 

Through the reflection coefficient a, the Kedem-Katchalsky formalism re- 
solves the competition between solvent and solute being transported through 
a shared pathway, e.g., a cotransporting channel permeable to both the so- 
lute and the solvent. The formalism itself, however, applies to any simple 
transport problem, regardless of whether a cotransporting channel is present 
or not (47J. In the latter case, not all three parameters {L p ,u,a} are inde- 



pendent. It can be shown 20] that in this case, a can be written as: 



(20) 



Introducing Eq. ( l2"Uj) into Eqs. ( TTTfT9"j) . one obtains the transport equations 
for the case where the solute and the solvent do not compete for the same 
cotransporting channel, e.g., in the case where they both diffuse through the 
phospholipid bilayer. Along the way, we will use the following simplifications: 
a Ac s ^> <p AII n w 0, Ap = 0, and introduce the notation more appropriate 
for describing the experimental setup: J v = (1/A) dV/dt, h s = (1/A) dN s /dt, 
P s = wa. Here, V is the total volume of the internal compartment, A is the 
area of the membrane, N s is the number of moles of permeating solute inside 
the internal compartment, and P s is the permeability of the membrane for 
the permeating solute. Using the described simplifications, the transport 
equations (j!7|18p can be written as 

^ = —L p aA [« + c e s ) - (4 + 4)] + P s V s A(c e s - c£) , (21) 
dA^ 

= -P s V s c s A(c e n - <) + P s A{c e s - c\) . (22) 
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Two terms contribute to the volume change in Eq. (1211) : the first one corre- 
sponds to the transport of solvent due to the osmotic pressure gradient, and 
the second one to the transport of solute. For dilute solutions, the second 
term is much smaller. Similarly, in Eq. (122]) . the first term is proportional to 
V s c s , which makes this term negligible for dilute solutions. 

As an example, we will apply the two-parameter model to the case which 
corresponds to the experimental setup: a single vesicle filled with a non- 
permeating solute is transferred to a reservoir filled with an isotonic solution 
of a permeating solute. Even though the solutions are isotonic, permeating 
solute diffuses into the vesicle, thus causing an osmotic pressure, which is 
in turn balanced by the inflow if water. Two cases can be distinguished: a 
flaccid vesicle changes its shape and becomes ever more spherical, while a 
spherical vesicle has already reached its limiting shape and can only grow 
by stretching the membrane. In this case, c e n = 0, c e s = const., and initially, 
c\ = 0, c l n = d nQ ; which is equal to c e s . 

In the case of a flaccid vesicle, membrane area A is constant, and the 
system can be characterized with the vesicle volume V and the amount of 
permeating solute in the vesicle interior, N s . Using the simplifications V w ~ 
V, which is appropriate for dilute solutions, Eq. ( 121]) and Eq. ( 122|) transform 
into: 

£ = -v*(<- (23) 

Here, V\ n \ denotes the initial volume of the vesicle, and c l n0 denotes the initial 
partial concentration of the non-permeating solute inside the vesicle. Often, 
the permeability of membrane for water is expressed as water filtration co- 
efficient Pf instead of hydraulic conductivity L p , the two quantities being 
bound by the relation Pf = aL p /V w . Alternatively, one can use the reduced 
volume v = V/Vq, instead of V, with Vq = A 3 ^ 2 / being the volume of 
a sphere with an area equal to A: 
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Here, v = V ini /V is the initial reduced volume of the vesicle. 

The other case is a spherical vesicle. In this case, neither its area A nor 
its volume V are constant; both can be, however, expressed in terms of the 
vesicle radius R, which is a convenient parameter in this case. Substituting 
A = AttR 2 and V = AnR 3 /3 into Eq. and Eq. ([21]), one obtains: 



dR - PV 

del 3P. 



Rq 
R 



dt R 



(27) 

(<5 ~ 4) • (28) 



Here, Rq is the initial vesicle radius; Vo = 47t.Rq/3. 

As figure |2] shows, both flaccid and spherical regimes can appear in the 
course of a single vesicle transfer. The two-parameter model allows us to 
mimic the same behaviour. Figure H] shows a time course of the radius of 
a vesicle cross-section upon transfer into a slightly hypertonic solution of a 
solute which can permeate the vesicle membrane. A spherical vesicle with 
Rq = 10 [im is initially filled with a 0.1 mol/L solution of a solute which 
cannot permeate the membrane, then transferred into a 0.105 mol/L solution 
of a solute which can permeate the membrane (P s = 2 • 10~ 8 m/s). Other 
parameters used were V w = 1.806 • 10~ 2 L/mol, Pf = 2.23 • 10~ 4 m/s. Upon 
transfer into a hypertonic solution, two processes are competing: the efflux 
of water, driven by the osmotic pressure mismatch, is the quicker of the 
two, while the influx of the solute, accompanied by the concomitant influx of 
water, is the slower one. Therefore, we first have to solve the system defined 
by Eq. flU} and Eq. fl23J, with V(0) = V Q and JV s (0) = as initial conditions. 
A deflated vesicle deforms into an approximate oblate spheroid; since in the 
experimental setup, the optical axis is usually aligned with the symmetry 
axis of the spheroid, the cross-section radius increases with the decreasing 
volume. Introducing v = V/Vq and x = R^/Rq, R\ being the spheroid visible 
semiaxis, one obtains the following relationship between v and x: 



v 



2 



1 + s/l -v 2 /x t 



2x 2 + ; In V ' = 4 . (29) 



, r i 



y/1 - V 2 /X 6 U - y/1 - V 2 /X S 



After a certain time, the vesicle becomes spherical again. From this point on- 
wards, it starts to grow while maintaining a spherical shape, and its behaviour 
is governed by Eq. ( 1271) and Eq. ( 1281 . with initial conditions R(0) = Rq and 
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Figure 4: The radius of a vesicle cross-section upon transfer into a hypertonic solution 
of a solute which can permeate through the vesicle membrane, as calculated by the two- 
parameter model. Calculation parameters are given in the text. 



c*(0) equal to the concentration of the permeable solute, which diffused into 
the vesicle during the previous step. 

3.3. One-parameter model 

Vesicle swelling is a two step process: in the first step, permeating solute 
diffuses into vesicle interior and thus causes osmotic non-equilibrium, and 
in the second step (which constitutes the bulk of swelling), it is followed 
by the water influx, which balances the osmotic non-equilibrium. Eq. ( |23|) 
and Eq. (124"]) can be regarded as a kinetic system. In the systems treated 
here, however, the influx of water is a much faster process than the solute 
diffusion, and thus Eq. ( |24"1) represents the rate-limiting step. In the case 
where solute diffusion is a much slower process than the influx of water, a 
simpler description can be obtained. 

In this simpler description, Eq. (|24"|) is retained, while Eq. fl23|) is replaced 
by an instantaneous osmotic equilibrium: 

<4 = ^ • (30) 
With c e s and N n being constant, this leads to dN s /dt = c e s dV/dt. Taking this 
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into account, and inserting Eq. (13T)|) into Eq. (I2"4"|) . one obtains 

(31) 



dV = P s AN n 
dt c5 V 



Eq. ( 131]) is general in the case that both A and V may vary with time. Again, 
two special cases can be considered: flaccid vesicles, where A = const., and 
spherical vesicles, where both A and V can be expressed in terms of the 
vesicle radius R. In the latter further simplification is possible by 

assuming that initially, the vesicle was spherical and in osmotic equilibrium, 
i.e., c e s = N n / (4nRl/3). In this case, Eq. (l3T|) can be rewritten into 

dR Ps(^)\ (32) 



dt \R, 
which can be readily integrated: 

/ 4PA 1/4 
R(t) =Ro\l + -^-) ■ (33) 

As long as deformations are small, R exhibits linear growth, R(t) « i?o + Pst- 
In this model, which assumes that the membrane permeability for water ex- 
ceeds the permeability for solute by such a large margin that water transport 
can be considered instantaneous, the only parameter governing the vesicle 
volume change is P s . Despite its simplifications, the one-parameter model 



has proved to be useful in certain situations |48|, |49 



3.4- The influence of membrane elasticity 

The treatment of osmotic swelling of spherical vesicles presented so far as- 
sumes that the membrane is infinitely "soft" and does not oppose its stretch- 
ing as the vesicle swells. In reality, the stretching modulus of a phospholipid 
membrane is finite. This means that the apparent value of membrane per- 
meability derived from the experiments with osmotic swelling of spherical 
vesicles is slightly underestimated: because the vesicle membrane opposes its 
stretching, the radius increases slightly less in a given interval of time than 
it would if the limiting factor was the membrane permeability alone. 

It is common to assume that the membrane area has a certain relaxed 
area Aq, and expand the free energy of membrane stretching around this 
value: 

W=^(A-A o y. (34) 
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Here, K is the membrane stretching modulus, and A = A^R^. We consider 
now the work p dV, needed to increase the radius of a spherical vesicle by a 
small amount dR. The corresponding change in the free energy is 

dW = K^-^dA . 

Substituting dV = 4iiR 2 dR, dA = 8nRdR, and equating pdV with the 
change in the free energy, one obtains: 

2K R 2 - Rl 

P=^^r Jl - (35) 

it JtCq 

In the case of osmotic equilibrium, this pressure is balancing the difference 
in the osmotic pressure: p = AIL 

When the pressure exerted by the membrane is taken into account, Eq. ([53]) 
transforms into 

^ = -L P A(AU - Ap) . (36) 

Here, An = a(c e — c*) and Ap = p e — p % = — p, where c e = c e s , c 1 = 
d n + c\, and p is defined by Eq. (135]) . In the system defined by Eq. ( ]27]) and 
Eq. (|28]) . the vesicle membrane would in theory grow indefinitely because 
the osmotic pressure is never entirely balanced. A membrane which opposes 
vesicle swelling limits the extent of swelling by the condition An = Ap. 

Rewriting Eq. ( |36l) for a spherical vesicle, one obtains an equation anal- 
ogous to Eq. (]28]): 



dR _ 

dF-" Lf 



a 




(37) 



For realistic parameter values, membrane stretching only adds a minor cor- 
rection to the vesicle swelling rate. Figure [5] shows the swelling of a spherical 
vesicle upon transfer into an isotonic medium of permeating solute (P s = 
2 • 10~ 8 m/s), obtained by the numerical solution of the two-parameter sys- 
tem defined by Eq. ([37D and Eq. ([28]), with R(0) = R and 4(0) = as 
initial conditions. The upper curve disregards the membrane stretching en- 
ergy (K = 0), while the lower uses the value K = 0.23 Nm _1 . Other 
parameters used in calculation are c e s = d n0 = 0.1 molL" 1 , R = 20 \un, 
L p = 1.645 • 10" 13 ms" 1 Pa" 1 , and a = RT = 2477 J mol" 1 at T = 298 K. 
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Figure 5: The radius of a vesicle cross-section upon transfer into an isotonic solution of a 
solute which can permeate vesicle membrane, as calculated by the two-parameter model. 
The lower curve shows the rate of vesicle growth when the finite membrane stretching 
modulus K is taken into account, while the upper curve shows the dependency for K = 0. 
Calculation parameters are given in the text. 



3.5. Swelling-burst cycle 

In reality, vesicles approximately follow the swelling course shown in Fig- 
ure [5] only until the critical strain of the membrane is reached. Once the 
critical strain is reached, the membrane ruptures, and the vesicle bursts. Sub- 
sequently, the membrane is resealed and another cycle of swelling commences. 
Experimentally, long trains of ~ 50 vesicle bursts have been observed. 

The phenomenon has been predicted and thoroughly worked out from 
the theoretical point of view by Kozlov and Markin 50| ; apparently unaware 



of their work, other authors attempted the same decades later [51] . Experi- 



mentally, the effect has been observed on erythrocytes when their suspending 



medium has been exchanged with a hypotonic one 52], on giant adhering 
vesicles (53[, on giant vesicles with equinatoxin-II-induced pores j^H, on giant 
vesicles with photoinduced pores 55[ , on giant vesicles with mellitin-induced 
pores 56|, and on giant vesicles upon transfer from a solution of a non- 
permeating solute to a solution of a permeating solute [37 . 

Kozlov and Markin [50 1 ] consider a vesicle filled with an osmotically active 
solute [i.e., non-permeating), and write an expression for the osmotic influx 
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J in f of water: 

J inf = —L P A(AU - Ap) , (38) 

where L p is the permeability of membrane for water (hydraulic conductivity), 
A is the membrane area, AIT = IF — IT* is the difference in osmotic pressure, 
and Ap = p e — p % is the difference in the hydraulic pressure. Eq. (138]) is 
identical to Eq. (J3E]1 above. 

The vesicle is treated as spherical, its radius denoted by R. Due to 
osmotic swelling, R > R , where Ro denotes the value of a non-expanded 
vesicle. Denoting R = Rq + AR, the corresponding changes in the vesicle 
volume and membrane AA/Aq « 2AR/R and AV/V « 3AR/R , 

respectively. A stretched membrane creates a pressure inside the vesicle - cf. 
Eq. (155|) . which opposes the swelling, 

2KAA 

Ap = , 

R A 

with K being the membrane stretching modulus. The elastic free energy 
associated with membrane stretching is given according to Eq. fl2! 



A pore present in the membrane lowers AA and thus reduces this energy 
term. However, on the other hand, creating a tension pore increases the 
total energy by the edge energy, equal to the product of the lenght of the 
pore edge, 27rr, and the linear tension 7: 

W p = 27r 7 r (39) 

Here, r refers to the radius of a circular tension pore. Eq. ( 139]) is an ap- 
proximate expression, valid for large pores. A consequence of a pore having 
a non-zero energy is that the intravesicular pressure does not drop to zero 
after the membrane ruptures. Instead, a residual intravesicular pressure re- 
mains: 

A Pr = . (40) 

This residual intravesiculare pressure causes the efflux of the intravesicular 
solution once the vesicle ruptures. The efflux J e s is approximately equal to 

J,« = ^ . (41) 
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Here 77 denotes the viscosity. 

Substituting Eq. ( l4"Uj) into Eq. (I3"oT) and integrating it, one obtains the 
time course of vesicle volume upon transfer into a hypotonic medium: 



Here, a = RT. The elastic energy of a vesicle with a pore is the sum of the 
membrane stretching energy and the linear energy of the pore: 



At a given AA, the area of the tension pore A p adapts in such a way as to 
minimize the total energy given by Eq. (T4"3"j) . If W is plotted against pore 
radius r for realistic values of K and 7 and for different values of AV/Vq, 
one can prove that for low values of AV/Vo, W is a monotonously increasing 
function of pore radius r (figure [6]). If AV/Vq is increased, we reach a critical 
value AV /Vq at which W(r) has an inflection point. By increasing AV/Vq 
even further, a local minimum becomes a global one, meaning that at some 
non-zero pore radius, the energy of the vesicle is lower than at r = 0, and 
by overcoming the energy barrier, the vesicle can jump from a poreless state 
into a state with a pore. 

Kozlov and Markin show that based on the difference between the influx 
and the efflux, given by Eq. fl3"5|) and Eq. (T4"T|) . three different scenarios are 
possible: 

Jiai > Jeff instant efflux of the whole vesicle interior 
Ji ri f = J c g steadily open pore 
Jinf < Jeff pulse-wise regime 

Kozlov and Markin conclude that for realistic parameter values, the pulse- 
wise regime is the most probable one, which is consistent with our own ob- 
servations. 

In the pulse-wise regime, it is possible to infer some relationships between 
relevant quantities during successive bursts by employing two assumptions: 

1. During the swelling phase, the amount of the non-permeating solute 
inside the vesicle remains constant. 

2. At burst, partial concentrations of both the permeating and the non- 
permeating solute remain constant. 




(42) 




(43) 
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Figure 6: The elastic energy of a vesicle with a pore plotted against reduced pore radius 
r/Ro for several different values of AV/Vq. Vesicle elastic energy is given by Eq. ([43]) . 
with i? = 5 urn, K = 240 mN/m, 7 = 20 pN. 

The first assumption is valid if the characteristic time for solute exchange is 
long compared with the time course of the experiment; the second assumption 
assumes that the interior of the vesicle is well mixed. 

While the validity of the above assumptions is more general, some simple 
formulas can be obtained for the linearized one-parameter case. Let us con- 
sider a vesicle in the pulse- wise regime (figure [7]). Due to osmotic swelling, 
the vesicle radius R increases from its relaxed value Rq until a critical value 
R c is reached, at which the critical strain for the membrane is reached, upon 
which membrane ruptures, ejects a part of the internal volume, and reseals 
in a relaxed state. As the (A c — A )/A amounts to a few percent, Eq. (13"3"1) 
can be linearized. 

Eq. ( 13U|) is valid only up to the first burst. Deriving it, we assumed 
c e s = N u /Vq. We can derive the concentration of non-permeating solute 
in subsequent swelling cycles using the assumptions above. Denoting the 
amount of non-permeating solute inside the vesicle after n bursts by Nn , 
the concentration of non-permeating solute inside the vesicle after n bursts 
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Figure 7: In the pulse-wise regime, vesicle radius R increases due to osmotic swelling from 
its relaxed value Rq until it reaches its critical value R Cl at which the critical strain for the 
membrane is reached. At that point, membrane ruptures, vesicle ejects part of its interior, 
upon which the membrane reseals and the radius returns to its relaxed value. 

by Ch \ Vq = 47t_Rq/3, and V c = 47ri?j?/3, one can write: 



A general expression for ch^ is then 



Vo 

Vc. 



(44) 



In our experimental setup, c£) = c e s holds. Substituting c l n with in the 
derivation of Eq. ( 1331) . one obtains the temporal dependence of vesicle radius 
after n bursts: 



R {n \t)=R Q 



1 + 



.1/4 



(45) 



Or, linearized, R {n) {t) = R + (V /V c ) n P a t. The time between two successive 
bursts is then n 

fVo\ -^c — Rq 

n ~\v c ) 

Denoting a c = (A c — A )/A and expressing V c /V w 1 + (3/2)a c , R c /Rq ~ 
1 + (l/2)a c , one obtains after linearization for the time elapsed between the 
n-th and the n + 1-th burst: 



Rq / a c 3a c 

t n ~ — H 1 n 

PA 2 2 



(46) 
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In this linear model, the time elapsed between subsequent bursts increases 
linearly with n. 

As the formation of a tension pore is a stochastic event, R c is not a 
constant, but can vary from one burst to another. To account for this, 
Eq. (TJ3]) can be written as 

*>=4?>"fe, (47) 

3=0 Vc 

(i) 

where the index j denotes each of the n bursts. We have introduced V c = 
AnRc /3, the volume of vesicle at which the critical strain for the membrane 
is reached in the j-th burst. 

3. 6. Critical phenomena in osmotic swelling 

Upon transfer of vesicles from a solution of a non-permeating solute into 
a solution of a permeating solute, vesicles often exhibit a transient increase 
of cross-section radius (figure |2]). We attribute this behaviour to a slight 
hypertonicity of the target solution, which causes the vesicles to slightly de- 
flate and change their shapes from an initially spherical shape into an oblate 
spheroid. The concentration gradient, however, drives the permeating solute 
into the vesicle interior, and the ensuing osmotic pressure difference serves 
as a driving force for water, which causes the vesicle to re-inflate. The infla- 
tion process has two stages: initially, the vesicle is flaccid and approaching 
the spherical shape, while in the second stage, the vesicle is spherical, its 
membrane being increasingly stretched. 

It has been shown that an inflating vesicle reaches the end of the first stage 
critically, through a continous transition. This was initially demonstrated for 



unstretchable membrane [57] . and later the theory has been extended in order 



to account for membrane stretching 43] . We shall provide a brief outline of 



the underlying theory, which employs a thermodynamic framework [58| . 

Let us consider a vesicle with Q entrapped molecules of a non-permeating 
solute. We assume that the vesicle is brought into a thermal equilibrium with 
the surroundings, characterized by temperature T and pressure po. Both the 
vesicle volume V and the surface area A are treated as free, independent 
thermodynamic variables. The Gibbs free energy of the system contains 
volume and surface contributions: 

G = G m + F 2D . (48) 
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Considering the enclosed solution as ideal and dilute, one can write 

G 3D = k B TQ (in ^ - + Po V . (49) 

The surface terms comprise the stretching and the bending parts, 

F 2D = F S + F b . (50) 

The stretching part is assumed to have a minimum at some value Aq] a 
quadratic expansion is used: 



P K(A-A 

r s — "7T 



,2 



2 A 

In regarding the bending free energy, only the contribution from undulation 
entropy is included, which describes the suppression of bending fluctuations 
as the vesicle shape approaches a sphere. Assuming small fluctuations around 



a spherical shape, it is given by [57 



N 

F b =--k B T\n(l-v). (51) 

Here, N is the total number of independent bending modes contributing to 
the membrane thermodynamics, and v is the vesicle reduced volume. It is 
important to note that through the reduced volume (which depends on both 
volume and surface area), F b provides the coupling between the volume part 
(t3d and the surface part i^D- The number of modes N is not known in 
advance and is a parameter extracted from the experiment. 

Using Eqs. |48H5T| the Gibbs free energy of a vesicle is defined via (T, p , Q,N), 
as well as A and V. By minimizing G with respect to A and V the equilib- 
rium free energy can be obtained. Introducing an intensive area a = A/N 
and its relaxed value ao = Aq/N along the way, one obtains a set of equations 
from which a and v can be calculated: 



2 1-v 

2 5n a(a — ao 

3 5 K a 2 





3/2 

) " = 


Q 


a 0/ 


Q c 


a x 


3/2 
) V = 


Q 



(53) 
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In Eq. ([52]) and Eq. ( 15"3"|) . we have denned 



g c = ?^, where V = ^1 , (54) 

and parameters 5jv and the first one related to the finite vesicle size, and 
the second one to the finite membrane stretchability: 

5 N = ^ociV" 1 / 2 , (55) 

5 K = ^LocK- 1 . (56) 

Two limiting cases of the described systen can be considered. Setting 
iV — > oo, 5n — > while keeping the stretching modulus finite brings us to 
the usual thermodynamic limit of an infinite system. In this limit, Eq. ( 152]) 
and Eq. (153]) are degenerate and yield the expected equilibration of the in- 
ternal and the external pressure: QksT jV = Po- The second limiting case 
corresponds to an unstretchable membrane: K — > oo, 5k — > 0. In this case, 
Eq. ( 153]) requires a = a , while Eq. ( 152]) yields the behaviour studied in 
Ref. & with v exhibiting a critical behaviour as N — > oo. 

Away from either limit, Eq. (|52|) and Eq. ( 153]) yield 

u 4 a(q - op) 
ffjc l^ = • (57) 

We can see from Eq. (|57p that attaining a perfectly spherical shape (v — 1) 
would reqire infinite strain a — )■ oo, which is not realistic. The deviation 
of the area from its relaxed state, a — a Q , is inversely proportional to the 
deviation of the vesicle shape from a sphere, 1 — v. As 5k "C 1, it takes a 
highly swollen vesicle, 1 — v ~ 5k, to produce a significant membrane strain 
(a - a )/a . 

In order to provide an adequate description of the transition, we define 
the control parameter q = Q/Q c — l, proportional to the number of enclosed 
molecules, and order parameter M = 1 — v , which serves as a measure of the 
deviation of vesicle shape from a sphere. Solving Eq. (152]) and Eq. ( 153]) for 
M in the vicinity of the transition, one obtains 




A ■ \/28 N + ~5 K . (58) 
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For q < 0, M(q) ~ \q\, which is appreciable, while for q > 0, M(q) ~ A/(4g), 
which is very small. The transition occurs over a region, determined by A. 
Thus, both the finite vesicle size and the membrane stretchability contribute 
to widening of the transition. In the limit A — > 0, i.e., for an infinitely large 
vesicle enclosed by a non-stretchable membrane, dM/dq has a discontinuity 
at (q — 0, M — 0). Setting 5 k = while keeping 5 at finite yields back the 
results obtained in [&7j | . Which of the two factors entering A is dominant, 
5 at or depends on their ratio 

5 N K/R 



<5k Po 

Using the values from our experiments - Rq ~ 20-50 (xm, Co ~ 0.1-0.2 mol/L 
- as well as K = 240 mN/m [59j, we obtain 5^/5k ~ 0.01-0.1. Therefore, 
in our conditions, the transition width is governed by the finite stretching 
modulus, and can be approximated by 

While Eq. ( 1581) fully describes the law of corresponding states for osmotic 
swelling of nearly spherical vesicles, we would like to transform it into a 
form which would allow us to compare it with the experiment, i.e., with the 
function R(t). 

The simpler part is the transformation of the time axis. As the concentra- 
tion of the permeating solute inside the vesicle is approximately two orders 
of magnitude lower than its concentration in the surrounding medium, the 
former can be neglected, yielding dQ/dt ~ PAqCq, which further leads to 

q = (3P/R )t + const. (60) 

Eq. (|60|) gives a simple linear relationship between the time t and the control 
parameter q. 

The transformation for the radius is less straightforward. Here we only 



reproduce the result, which is worked out in detail in (43|. Denoting g exp = 
(Ri(t) / R ) 3 — 1, where Ri(t) is the radius of vesicle projected cross-section, 
one arrives at the following scale and shift transformation, 



/(0 



15A 



icxp 




eA + ^(l5A 2 ) 1/3 



(61) 
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where £ = q/A and q = (3P/Ro)t. Using the permeability P and the transi- 
tion width A as two fitting parameters, one can verify that the experimental 
data for all recorded vesicle transfers collapse onto a single universal function, 

/(£)= (v / l + e-e) V2 • (62) 

Fitting Eq. ([6"2"j) to the recorded data R\(t) of vesicle transfer into an isotonic 
solution of urea yields the permeability of membrane for urea, P s = 0.013 ± 
0.001 (xm/s. The obtained permeabilities for glycerol and ethylene glycol, on 
the other hand, exhibit concentration dependence. The values obtained for 
glycerol range from 0.0053 at c = 0.1 M to 0.019 ± 0.006 [im/s at 0.2 M, 
and those for ethylene glycol range from 0.046 ± 0.006 |j.m/s at 0.1 M to 
0.085 ± 0.01 nm/s at 0.2 M. 



4. Experimental analysis 

In this section, we will present a method for determining membrane per- 
meability for a given solute from the analysis of swelling-burst cycles [39( , and 
demonstrate that its results are consistent with the analysis of the osmotic 



swelling of flaccid vesicles |43 



In section [31 Eq. (145]) has been obtained, which, when linearized, shows 
that before the first burst, the slope of R(t) equals to the permeability P s , 
and after n bursts, the slope becomes {V /V c ) n P s if all bursts are equal, and 

Ili=Q (^o/^c ) if they are not. The underlying assumptions are that the 
vesicle is spherical at the onset of the first swelling cycle, and the concentra- 
tion of the non-permeating solute in its interior is equal to the concentration 
of the permeating solute outside. 

These assumptions are not always met. In 37 out of 47 recorded vesicle 
transfers in our experiments, the dependence of the vesicle cross-section ra- 
dius as a function of time exhibits an initial transient maximum or "bulge" 
like the one seen in Figure [2j In those cases, the vesicle is filled with a mix- 
ture of a non-permeating and permeating solute already at the onset of the 
first swelling cycle, and Eq. (|45p needs a small correction to take that into 
account. 

An estimate of the mole fraction of the non-permeating and the perme- 
ating solute at the onset of the first swelling cycle can be obtained from the 
initial transient maximum. If we approximate the shape of a partially de- 
flated vesicle with an oblate spheroid and denote x = Ri/Rq, where R\ is 
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the radius of the vesicle cross section at the crest of the "bulge" and R is 
its value in the "valley", i.e., at the onset of the first swelling cycle, we can 
compute the reduced volume of the vesicle (v) using Eq. fl29|) . 

We can show that the reduced volume of the vesicle is equal to the mole 
fraction of the osmotically active solute at the onset of the first swelling cycle. 
At the crest of the bulge, the vesicle is flaccid with a volume V = vV < V , 
and containing iV n moles of non-permeating solute. At the onset of the first 
swelling cycle, the vesicle is spherical with a radius Ro, and containing N n 
moles of non-permeating solute and N s moles of permeating solute. As the 
rapid osmotic exchange of water ensures that the total concentration is equal 
in both cases, NJV = (N n + N B )/V , it follows 



x n 



N n + N s 

We can amend Eq. (1451) with this correction, yielding 



R (n \t)=R 



VoV 4x n ,P^ 1/4 



V, Rn 



(63) 



We need to note that taking the value of x n at the onset of a swelling cycle 
is an approximation, as x n actually changes during the course of a swelling 
cycle. 

Treating the vesicle swelling-burst cycles as piecewise linear, the slope of 
the fc-th cycle can be expressed as 

ARM ^^tJ. Vc 



3=0 Vc 

While membrane permeability for the permeating solute can be in principle 
computed from an individual vesicle-swelling phase between two bursts, a 
more reliable estimate is obtained by averaging it over all n bursts in the 
cycle: 

1 x " r 1 Ai?^) k . l V 
fc=l n j=0 Vc 

Eq. (1651) gives the estimate for the membrane permeability based on a single 
recording of a vesicle transfer. Averaging over several recordings of vesicle 



F * k) II 77(7) ' ( 64 ) 
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Table 1: Permeability estimates for glycerol, urea, and ethylene glycol, determined from 
series of micrographs of giant unilamellar vesicles upon transfer from a sucrose/glucose 
solution into an isomolar solution of a given solute. "Burst train" refers to the analysis 
presented here, and "critical swelling" to an alternative analysis on the same data set. 





Permeability P s [|xm/s] 


solute 


"burst train" 


"critical swelling" 




0.1 M 0.2 M 


0.1 M 0.2 M 


glycerol 


0.0077±0.0009 0.016±0.003 


0.0053 0.019±0.06 


urea 


0.014 ±0.001 0.013±0.001 


0.013 ±0.001 


ethylene glycol 


0.054 ±0.005 0.10 ±0.01 


0.046±0.006 0.085±0.01 



tranfer are needed to obtain a more reliable estimate of membrane perme- 
ability. 

Table [1] summarizes the permeability data for glycerol, urea, and ethylene 
glycol, obtained by using two different analyses. The term "burst train" 
refers to the analysis presented here, while the term "critical swel ling " refers 
to the alternative analysis presented in the preceding section [43], where 
vesicle transition from a flaccid to a spherical state before the first burst is 
observed. 

One can see that in the case of urea, the estimate for membrane per- 
meability is independent of solute concentration, which is consistent with 
our finding in 43j. On the other hand, both polyols exhibit a concentra- 



tion dependence of the (apparent) membrane permeability, which is again 



in agreement with our own recent findings [43||. The trend is even more ap 



parent in Figure [SJ A possible explanation for the concentration dependence 
of the apparent membrane permeability is strong evidence for the affinity 



of polyols to phospholipid headgroups [60|, |61|, which results in the solute 
concentration adjacent to the membrane being larger than in the bulk solu- 
tion, which consequently leads to an apparent permeability coefficient which 
is larger than the actual one. This affinity of polyols to phospholipid head- 
groups may also account for a large scatter of published data on membrane 
permeabilities for polyols. 
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Figure 8: The apparent permeability of POPC membrane for glycerol, urea, and ethylene 
glycol, obtained by the "burst train" analysis. Fig. 6 in Ref. [43j shows analogous results 
obtained by the "critical swelling" analysis. 
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5. Discussion 



The existance of reproducible techniques for producing giant unilamellar 



vesicles (GUVs) [35|, |62|, the availability of phase-contrast microscopy and 
high-resolution digital cameras, along with the efficient procedures for au- 
tomated extraction of vesicle contour from the images and determination of 



the radius of vesicle cross-section |39j, render the determination of membrane 
permeability of GUVs a viable alternative to the established techniques em- 
ploying either dynamic light scattering on sub-micrometer vesicles or planar 
lipid membranes. 



GUVs were used before in studies of membrane permeability [32|. The 
study by Boroske et ai, however, was designed to observe osmotic shrinking 
of GUVs, while in our study vesicles swell osmotically (with the exception of 
the initial phase, where a transient osmotic shrinking may be observed). Os- 
motic swelling of vesicles in a controlled manner guarantees a simpler spher- 
ical geometry, through which both the vesicle area and its volume are ex- 
perimentally accessible variables. By contrast, in the experiment of Boroske 
et ai, vesicle volume was the only experimentally observable variable, and 
the authors had to come up with a plausible explanation for the apparently 
"missing" membrane area. 

On a very superficial level, the continuous transition between the "iron- 
ing" and the stretching regime of a vesicle which swells osmotically due to 
the permeation of a permeable solute might seem similar to the "minimal 
volume" technique employed by Sha'afi et ai (cf. Fig. 3 in [19| and Fig. 1 



m 



43)). However, one needs to be aware of the fact that Fig. 3 in 19[ shows 



vesicle volume as a function of time, while Fig. 1 in |43j shows the radius 
of a non-spherical vesicle cross-section, and throughout the course shown in 
the diagram, vesicle volume is monotonously increasing. In our experiment, 
minimal vesicle volume appears at the top of the initial bulge (visible around 
40 s on the graph in figure [2]) . We have not attempted an analysis analogous 



to the one performed by Sha'afi et al. [19[ or Hill and Cohen [22| . because the 



initial moment in which the vesicles are brought in contact with the target 
solution is not defined as well as in the stopped flow experiments. The au- 
thors of earlier studies, on the other hand, seem to be unaware of the critical 
transition between the "ironing" and the stretching regime in the course of 
osmotical swelling of a vesicle. 

While offering a direct visualization of osmotic swelling, the method de- 
scribed here has its limitations. The three solutes tested here - glycerol, 
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urea, and ethylene glycol - induce significant changes in vesicle size on the 
timescale of minutes. With solutes which are much more or much less per- 
meable, the rate of changes may be either too fast for video camera, or so 
slow that other effects start to dominate the experiment (e.g., convection 
in the measuring chamber). Another limitation is imposed by the refractive 
index of the solution. The refractive index of the target solution in general 
differs from that of the vesicle interior, and the difference between the two 
refractive indices increases with the increasing concentration of solutions. 
At high enough concentration, the halo around the vesicle becomes blurred, 



and the algorithm used to determine vesicle contour [39] fails. With our 
present equipment, we estimate that concentrations significantly higher than 
0.2 mol/L are outside our reach. 

Finally, we want to compare the results of two techniques dubbed "burst 
train" and "critical swelling" in Table [TJ As can be seen, they produce con- 
sistent results which are within the experimental error of each other. The 
experimental error estimates for both methods are also approximately the 
same. The "burst train" technique is more forgiving towards minor glitches, 
such as vesicle temporarily drifting slighly out of focus, which results in an 
overestimate of vesicle radius (seen around 370 s in figure [2]). However, one 
needs to keep in mind that the "burst train" technique requires approxi- 
mately 5 times as much data as the "critical swelling", and consequently, 
minor glitches are more likely to occur during a longer period. We therefore 
lean towards recommending the "critical swelling" technique, which is less 
resource-greedy while giving the results with the same precision. Our rec- 
ommended procedure is therefore to keep recording the vesicle transfer up 
to the first burst (which is easy to notice during the experiment), and carry 



out the analysis desribed in Ref. |43 
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